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Abstract 

A calculational approach in fluid turbulence is presented. A tur- 
bulent environment is thought to be a condition of the system for 
which the time and space evolution of the dynamical variables are 
chaotic, but for which a statistical distribution function is applicable. 
The time-evolution of an initially given probability distribution func- 
tion for the dynamical variables is given by the deterministic evolution 
equations that govern the system, called deterministic chaos. A fluid 
environment is a dynamical system, having energy input mechanisms 
at long distance scales occurring on the spatial boundaries and energy 
dissipation mechanisms due to viscosity occurring on small distance 
scales. Hence, the fluid has an attractor, called the Strange Attractor 
in the turbulent regime, for its dynamical variables. It is proposed here 
that the fluid probability density functional also has an attractor for 
its time-evolution, and an approach to generating this time-evolution 
of the statistics is presented. The same mechanism that causes the 
dynamical variables to have an attractor in phase space, that is the 
tendency for the equilibration of energy input rates and energy out- 
put rates to set in, also causes an arbitrary initial statistics to evolve 
toward an attracting statistics, which is stationary in time. It is this 
stationary statistics that allow the Kolmogorov scaling ideas to have 
applicability. The evolution of the fluid's statistics can be set up as 
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part of an approximate space-time path integral. Ensemble averages 
of any dynamical variable can be formulated in terms of this path in- 
tegral. Fluid space-time configuration sampling techniques naturally 
suggest a useful way, using a relatively arbitrary initial statistics func- 
tional to calculate averages. 



1 Description of the Mathematical Ap- 
proach to Fluid Turbulence 

Let us set up the evolution equations for incompressible fluid dynamics. 
The extension of the proposal to be described to compressible fluid 
dynamics will pose no unusual difficulty 
We have, 

dv -> « _ 

p ~dl = +v V} 

where / is an external force density, such as due to a scalar pressure 
field, and rj is the coefficient of viscosity, v is the fluid velocity field. 
We also have, 

v. W + g = o. 

Here, p is the mass density of the fluid. If this mass density is a 
constant, then the velocity field is divergenceless. Also, 

dv dv _^ 



So we have the fluid dynamic system, 
dv VP 



Vv + vV 2 v (I) 



at P 

V-v = 0, (2) 

where P is the pressure, and v = 2 

We drop the external force density in what follows. This is not 
an essential step. What are needed are a set of interacting fields, not 
restricted to being velocity fields, together with stationary boundary 
conditions to allow the deterministic time-evolution of the set. We also 
associate with the velocity field a probability density functional, p[v, t]. 
The fluid statistics time-evolves according to deterministic chaos [lj [2j, 



P[ v f>tf]= / d[vo]K[vf,tf,v ,t ]p[vo,to], 
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where the kernel is, 

K[v,t;v ,t ] = S[v f - v dassica i[t f ;vo,t ]}. 

That is, the number, p, associated with the velocity field vo at time to 
will be associated with vj at time tf, where Vf is the velocity field vo 
deterministically evolves into from time to- Given a functional of the 
spatial velocity field, A[v], its ensemble average, at time tf is, 

<A[v] >= J d[vf]A[v f ]K[v f ,t f ;v ,t ]p[v ,t ]d[v }. 

We want to propagate the fluid's statistics according to determin- 
istic chaos, even though the detailed fluid orbits are chaotic. Let, 

<A[v}> = J A[v f ]p[v f ,tf]d[vf] 

= J A[ Vf ]K[v f ,v ]p[voMd[v f ]d[v ] 

= J A[f[v ]]p[v ,t ]d[vo], 

where, 

K[vf,v ] = 5[vf - f[v ]] 

= J S[v f - /i[vi]](5[ui - /i[u ]]d[wi] 

= J 5[vf - f2[vi]]S[v 1 - f2[v2]]S[v 2 - f2[vo]]d[v 1 ]d[v 2 ] 

= J 5[v f - f3[vi]]S[vi - f3[v 2 ]]5[v2 - f3[v3]]S[v 3 - f3[vo}}d[vi}d[v 2 ]d[v 3 } 

. . . , 

where the velocity fields, v±, v 2 , V3, etc. occur in chronological order, 
v\ being closest to time tf. Eventually, we have an Jm, where M is 
large, such that vm = /m[^o] is infinitesimally different from vq. 
So, 

<A[v] >= J d[v]A[v f ]S[v - f M [v]]p[v ,t ], 

where, 

d[v] = d[v f ]d[vi] ■ ■ ■ d[v M ]d[v ], 

and, 

< A[v] >= I A[f M [- ■ ■ [f M [vo]] ■ ■ -]p[vo,to]dM- 
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We must have, 



K[v f ,v ] = 6[v f - f[v }] = Slf-^Vf] - v ]. 

We have, 

p[vf,tf] = J S[f' l [v f ] - v ]p[vo,t ]d[v ] 
= p[v ,t ]. 

The exact rule, /a/M, requires a solution of the fluid dynamic equa- 
tions, incorporating boundary conditions. The exact rule, /a/M, is 
difficult to find. Let us use an approximate rule, Fm[v]- 
Then we may say, with motivation to follow, 

<A[v]> =Jd[v}A[v f }5[v-F M [v}}X[V-v}X[v-VB]p[v ,t }. (3) 

The functional integration is over all space-time velocity fields within 
the spatial system, between times to and tf. 5[v — Fm[v]] is a space-time 
delta functional. Fm H generates v from a v an infinitesimal instant 
earlier. The A functionals are evaluated at a particular instant. The 
needed properties of the A functionals are A[0] = 1, and X[g ^ 0] = 0. 
X[g] could be, 

lim e"f 2/ei . 

ei->0+ 

We have, then, 

< A[v] > = J d[vo]A[F M [F M [- ■ ■ [F M [vo]) ■ ■ -]p[vo, t ] 

X[F M [- ■ ■ [Fm[vq\] ■ ■ •] - v B ] ■ ■ ■ X[F M [v ] - v B ]X[v - v B ] 
A[V • F M [ ■ ■ [F M [vo]} •••]]••• A[V • FmM]A[V • v ). (4) 

The right-hand side of (@J equals < A[v] >, because p^o^o] will be 
non-zero only for i»o's that make all the A's equal one. This means that, 
ideally, we utilize only spatial velocity fields that are divergenceless, 
and that satisfy the stationary boundary conditions. 

The velocity field has an attractor, determined by the stationary 
boundary conditions on the fluid. When the boundary conditions al- 
low steady laminar flow to become established, the attractor consists 
of a single velocity field. When the Reynolds number becomes large 
enough, bifurcations set in, or the onset of instability occurs, and the 
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attractor begins to consist of more than one spatial velocity field. In 
the turbulent regime, the attractor consists of many velocity fields, 
and the fluid accesses these velocity fields according to a probabilty 
distribution. 

Given a functional of the spatial velocity field, A[v], and the fluid 
dynamic system of equations (J2J, we will say that its ensemble 
average when the system has reached its attractor is, 



< A[v] > = lim d[v]A[vf]5[v - F M [v]]\[V ■ v]\[v - vb] 

■p[vo,t ]. (5) 

The delta functional condition, 5[v — Fm[v]] implements equation 
Fm[v] is an approximate rule for carrying v from an earlier instant to 
a later instant. In the first approximation, 

F M [v] = v - At(v-Vv - iA7 2 t>), 

where v is evaluated at an instant At earlier. A[V • v] implements a 
zero divergence condition on the spatial velocity fields, and X[v — vb\ 
requires the spatial velocity fields to have values «b on the spatial 
boundaries. 

Let us consider the path integral © to be on a space-time lattice. 
We could use 

_<<£ 

8{x) = lim 



+0+ y/We 

We have for the average of A[v] in the steady-state (attractor), letting 
ei(Ax) 2 = e, where e\ is in the A functional for the zero divergence 
condition, 

< AM > = lim lim / d[v]e- H[v]/e A[v f }p[v }. (6) 

e-»0+ tf—to— >oo J 

Also, 

p[v ] = p[v ,t ], 

and the space-time integration measure, d[v], is with boundary effects 
neglected, 

d[v} = (-}=f N l[dv ljM . 

H[v] is a functional of the lattice space-time velocity field. Also, ne- 
glecting boundary effects, 

H[v] = ^2({vi-vi-i-g[vi-i]At) 2 +(v x ,ijki-v Xt i-ij k i+- • ■) 2 )+^{vijki-VB,ijki) 2 , 
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g[vi] = -v- Vvi + vV 2 v t . 



Or, 



r i (Vijkl — Vi-i jkl) {Vijkl — Vi-l jkl — Vi-l jM + Vi-2,jkl) . 

9 W = -v x , ml — +• • ■+„ ^ +• 

N is the number of space-time lattice points in the system. We have 
as a sum over all space-time lattice points in the system, neglecting 
boundary effects, and J2 as a sum over an space-time lattice points 
on the spatial boundary. Also, ijk are spatial lattice indicies, and I is 
the index in the time direction. 

This discretization technique is expected to get better as one in- 
creases the lattice fineness and makes use of higher order finite dif- 
ference approximations applied to partial differential equations A 
good approximation to the attracting statistics as a starting point will 
shorten the evolution time required for accurate results for averages 
occurring in the steady state. The path integral {JjjJ) can be evaluated 
with Monte Carlo techniques utilizing importance sampling. A calcu- 
lation of the stationary velocity field that would exist, for the given 
boundary conditions, if that field were stable, should be a good starting 
point from which to begin a sampling of the velocity field space-time 
configuration space. 
We have, 



< A[v] >= J 



~! n 

i=i 



j= e - H W e d[v], 



where, 



and the space-time configurations are distributed according to the 
weighting e~ H *- v '< € . n is the number of space-time configurations in 
the ensemble. -A[u/t] is the observable A[v] evaluated at the final time 
slice of the i th space-time configuration. The value p[i?o,i] is attached 
to the initial time slice of the i th configuration. 
We need, 

1 = p { Vo d 

^ n 

i=i 

So, we must do the integral, J, and constrain Pi v o,i] to equal n/J. 

Going to higher order finite difference approximations for the deriva- 
tives will eventually require consideration of space-time boundary ef- 
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fects for each space-time point in the system. Even for the initial ap- 
proximation, boundary considerations must be made for points near 
the boundaries. 



2 Summary 

We have said that the time evolution of the statistics also has an 
attractor for the following reasons; 

1. It is a way to get a solution to the problem of arriving at the 
steady-state turbulence statistics. One knows that the steady 
state statistics is stationary with respect to its time-evolver. Propos- 
ing that this statistics is the attractor of its time-evolver means 
one does not have to have the statistics to get the statistics, 
thereby offering an approach to the closure problem for the de- 
termination of correlations. 

2. The statistical physics approach has been successful in equilib- 
rium thermodynamics where the derivation of the microcanonical 
ensemble can be taken as the indication that equilibrium is the 
attractor for the dynamical system when the boundary condi- 
tions on the system input no energy. In the attractor, the mean 
energy input is equilibrated with the mean energy output, be- 
cause in the attractor dissipative losses have effectively shut off, 
and the system becomes effectively Hamiltonian. The stationar- 
ity of the statistics requires the vanishing of the Poisson bracket 
of the statistics with this Hamiltonian resulting in the statistics 
of equal a priori probability. 

3. The dynamical system, of which a fluid is an example, has an 
attractor 0j. The dynamics of the statistical approach should 
mirror the dynamics of the actual dynamical system. 

4. The statistics of the dynamical system prior to reaching the at- 
tractor has no reason to be unique. The statistics of the attractor 
is unique, in which the geometry of the system, the stationary 
boundary conditions, and the viscosities, all of which determine 
the Reynolds number, play a crucial role in determining the at- 
tractor. 

5. The stationary statistics of the fluid occurs when the equilibra- 
tion of energy input and energy output has set in p]. 
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3 Conclusions 



In the discretized version of the path integral that attempts to arrive 
at the stationary statistical effects in the generation of the ensem- 
ble average of a dynamical variable, using an approximate rule for 
the dynamics, one should arrive at, in the continuum limit, a greater 
insensitivity to the initial statistics and a generation of steady-state 
statistical effects. One is trying to use the turbulent transience to get 
at steady-state turbulence effects. These steady-state statistical effects 
become the backdrop for Kolmogorov's ideas of self- similarity and the 
resulting scaling exponents. 
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